Splitting Functions in the Light Front Hamiltonian Formalism 

Raj en Kundu* 
Saha Institute of Nuclear Physics, 
1 /AF Bidhannagar, Calcutta- 700 064, India 
(17 May 1995) 

In this paper, we have calculated the lowest order splitting functions occurring in the Altarelli- 
Parisi equation in the light front hamiltonian formalism. Two component perturbative LFQCD in 
light cone gauge is our guideline and we have used LFTD-approximation for dressed quark and gluon 
states. 



I. INTRODUCTION 

It is believed that QCD is the potent theory for the correct description of the strong interaction, but explaining 
the real situation like the hadrons is still far from being complete. In the low energy non-perturbative sector, it is 
plagued by the confinement and the vacuum structure problems among others, which must be addressed properly to 
have the true picture of the QCD bound state. Recently, there has been an increasing amount of interest to explore 
the non-perturbative domain using the LFQCD in the light cone gauge Q. Asymptotic freedom and the factorization 
theorems which separates the soft and hard part of a QCD process, enables one to use the perturbative treatment 
in the high energy sector. The soft part measures the low energy (nonperturbative) properties of quarks and gluons 
in the parent hadron and is connected to the parton distribution functions. The hard part is relevant for the scale 
evolution of the hadronic structure functions and we shall be mainly concerned with that in this paper. 

It is a well known fact that the light front current commutators have the bilocal structure in it and the fourier 
transform of such bilocal matrix element gives the structure function of the deep inelastic scattering It is just one 
step further to see that parton pictures emerge directly from it. Scaling of the structure functions is very evident in 
this calculation. The systematic violation of the scaling comes only when the QCD corrections are taken into account, 
as can be seen, for example, from the Alteralli-Parisi equation Q. 

Parton model is best realised in the infinite momentum frame. Any kind of calculation involving parton model 
in the equal-time framework assumes some infinite momentum limit (which is sometimes conceptually difficult) and 
imposes some physical constraint on the gauge field to get the meaningful result. We prefer to attack such problems 
using light front framework and light cone gauge. Thus, with all the arsenal of two-component perturbative LFQCD 
in the light cone gauge we study the problem from the first principle without any further assumption. In this paper 
we calculate the lowest order splitting functions (P qq ,PG q ,PqG and Pgg) ^ occurring in the AP-equation following 
this path. This illustrates the efficiency and both conceptual and calculational simplicity of perturbative LFQCD to 
produce some physical results. 

In this calculation we have used the light front Tamm-Dancoff approximation (LFTD) ||, to describe a dressed 
particle. This formalism directly gives the probabilistic interpretation of the splitting functions. Similar calculations 
of the splitting function have been done in the equal time formalism taking infinite momentum limit in the axial gauge 

i- 

In this paper, we have briefly reviewed the necessary results in the first part and we use those results in our 
calculation in the latter part of the paper. For detailed discussion one should consult the references mentioned within 
the text. 



II. A BRIEF REVIEW 
A. LFQCD hamiltonian 

In our calculation we use two component formalism of LFQCD. Here we write the LFQCD hamiltonian as a free 
term plus the interaction term: 
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H = J dx-d 2 x ± (H + H int ) (2.1) 



where 



'Hint — 7~Lqqg + H-ggg + H-qqqq + H-gggg (2-2) 

Here TLint is divided into four parts depending on the nature of the interaction which can be read off from the suffix, 
for example, H. qqg implies the interaction of two quarks with a gluon and so on. First two terms of the interaction 
are of the order of g, the coupling constant and the order of the other two is of g 2 . Since, in our calculation we shall 
only be concerned with the lowest order in the coupling constant, we present here the first two terms explicitly. 

= {-2^) {B±.A ± ) + (a.A ± ) (J^j (a.d ± +m)+(J^j {a.8 ± - m)(a.A ± ) J £ (2.3) 

Uggg = Q f*° {(^JA^ + ) (±) (2-4) 

Here .A* (x)'s and £(x) stands for the dynamical component of the gauge fields and spinor fields. They are the 
solutions of the corresponding equation of motion and given by the following expressions, 

Ai w = E / 2 tl*C] [eW "' x)e ~ l9X + h ' c] (2 ' 5) 

£(*) = £xa / [b{p,X)e-^ + d\p,-X)e^] (2.6) 

2 2 _i_ m 2 

with 9" = fr and P~ = Pj [p+7 • 
Here A is defined to be the helicity, 

^ _ | 1 for gluons and ^ _ | 5 for quarks. ^ 7^ 

The gluon polarization vectors are e\ = (1, i) and e % _ 1 = ^(1, — i) and quark spinors are x\ — ^ q 



Also the creation and annihilation operators follow the basic commutation relations 



[a(q, X),a\q', A')] = 2(2tt) V<*AV<5 3 (<7 - </) 

{b(p, A), 6t(p', A')} = {d{p, A), rft(p', A')} = 2(27r) 3 <5 A , A ^ 3 (p _ p >) (2. 
In our calculations we also require the multiple principal value prescription, given as 

(jf) ^ X '^ = (i) / dx l dx 2-- dx n e ( X ~ - x l)~~e( x n-l ~ X n)f( X n) 

' 1 1 + ^— ] f(k + ) 



2 \k+ +ie k+ - /< 
1 



[k+\ n 



f(k + ) (2.9) 



B. Two component perturbative LFQCD 

Here we present the rules that we have used to calculate the necessary matrix elements, 
a) Draw all topologically distinct x+-ordered diagrams. 
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b) For each vertex, include a factor of 16n 3 S(pf — pi) and a simple matrix element given below. Each gluon line 
connected to the vertex contributes a factor from the normalization of the single gluon state. 
Necessary matrix elements are as follows. 



H, 



ggg 



t / 2kl ( a -P2 ~ im ) i i ( a -Pi + im ) 

-gTi a x{r i { Pl ,p 2 ,k)xxA 



XA 2 e A 



(2.10) 



H, 



ggg 



l 9.J e Ai e A 2 e A 3 1 



(fc 2 - hy 



[kt\ 



+ 



(h + hy-^-ikt + k^Su 

Ml 



[k~2 k~z) 



-(k' + k 2 ) 1 



Hi 



(k+ + k+) 



Si i 



(2-11) 



Using energy-momentum conserving <5-function and expressing them in terms of relative momentum defined as x — 



k + 



k+ 

j^-fjPi, we can write 

r( P , P -k,k) 



= 2 



a- 7 (jP — ti) — im ; a aipi + im 
-<j — a 



[p+ - k+] 



r«'(p,A,p-fe) 



, f-. -K l - a z (a.Kj_) + ia z mx 



(p-2ky--^ T - ] (p+-2k+) 



(2.12) 



+ 



(fc - 2 P y 



[k+] 



{k+ - 2p+) 



K 3 Su + 



(2.13) 



C. LFTD approximation in (3+1) dimension 

Here we briefly review the light front Tamm-Dancoff (LFTD) approximation in (3+1) dimensions which will be 
necessary for our purpose. LFTD approximation, in its simplest form, is just the Tamm-Dancoff approximation 
applied to light front field theory. In Tamm-Dancoff approximation we usually truncate the fock space as much as 
necessary and plausible, and then expand the dressed particle state in the truncated fock space. For example, we can 
write the state corresponding to a quark as follows: 

\quark) = \ip{p + ,p ± )) 

- Clip) -^(p,x)\o) 



Ai A 2 



= ci^i + c 2 ?/>2 (say) (2-14) 

[Here A's stand for spin as well as the colour indices.] 

Thus we have truncated the fock space to a bare quark plus a quark and a gluon state, so that the maximum 
occupation number is two. For higher approximation we can take a three particle state also. Here ci and C2 arc 
the probability amplitude of finding a bare quark with all the momentum and of finding a non-interacting pair of 
quark-gluon with given momentum sharing. 
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Similarly, we can write a gluon as, 
\gluon) = \tp(p + ,p- L )} 

= / l(p) qt(p,A)|0> 



+ ]T I I ^^SHp-q-k)bHq,\ 1 )d!(k,\ 2 )\Q)c 2 ( q ,k) 



6 3 ( P - q - k)a\q, \ x )a\k, A 2 )|0)4(g, k)) 



dq + d 2 q ± f dk+d 2 k A 



^2(2^)35+ 7 v /2(2 7 r) 3 fc+ 



C\\bare — gluon) + C2\quark — antiquark) + c' 2 \gluon — gluon). 



(2.15) 



We also have to use within this fock space some momentum cutoff and in general we should have a step function 
such that the invariant mass be less than the cutoff. 

One last comment about the normalization of the state. It is to be noted that the state is not yet normalized. 
Using the normalization condition one can get rid of one of the c's, for example, ci(p). 



III. SPLITTING FUNCTIONS 



A. Calculation of P qq and Pc q 



First consider the state of a dressed quark according to the Tamm-Dancoff approximation. 

ci{p) hUn m) + sr f g^jV f dk+d^ 1 - 



b\q,X 1 )aHk,X 2 )\0)c 2 (q,k) 
= ci^i + c 2 ^ 2 (say) 

Here A's stand for spin as well as the colour indices. 



S 3 {p-q-k) 



Now, we consider the dispersion relation (P ~ + P~) = — p + ± , and operate this on the quark state. 



(p - + p-)\tP) 



p+ 
M 2 + P\ 
P+ 



(3.1) 



(3.2) 



Now take the projection of this equation onto non-interacting one-fermion and one-fermion-one-boson state. This 
leads to 

A,f2 , p2 

(p'\(Po+p v -M) = (p'\ ± W 



{p'k'\(p - + p-M) = {p'k'\ 



P^ 
M 2 +P 2 
P+~ 



(3.3) 



where, 



\P' 



1 



P >,k>)= 1 



^2(2^)3 v /2(2 7 r) 3 fc+' 
Now, from equation (|3.3|), we get, up to the lowest order in g, 



tf(p',\[)at(k',\')\0) 



(3.4) 



c 2 {p k \P \ip 2 ) +ci(pk \P V \ipi) = c 2 {pk 



P^ 
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or, c 2 
or, c 2 



M 2 -i- P 2 
(P'k'\ p + H ^ 2 )-{p'k'\PM 

M 2 + (p' ± + k' ± ) 2 



(p'+ + k'+) 



i-F +Pl 



k'+ 



d^'fc'IP-^). 

cx^'k'lP-^x). 



(3.5) 
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Now, (p'k'\P v can easily be calculated using LFQCD perturbation theory. We can use the rule already mentioned 
and take M = mp — vtlq = limit of the above expression. Also, take the total transverse momentum to be zero. 



=4> p' ± + k' ± = or, p' ± = —k' ± = k± {say) 
Thus, we parametrize the momentum as follows: 

total mom. p = (p,p, 0j_) 



(3.6) 



quark mom. p' = I (1 — x)p + 



k , 



2(1 -x)p 



gluon mom. k! = ( xp H —,xp, —k± 

' ' 2xp 



(1 - x)p, k± 



(3.7) 



With the above prescription, we can write equation (3.5) as, 



C-2 



1 



1 



x 1 — X 

^ , 2 (ll)V ^ C ' l(p/fc/|Arl ^ )|2 -2 ( 2^ + 



(3.8) 



[Note that in the last step, we have written normalization factors and the factor coming out of energy momentum 
conserving S- function separately, so that the matrix element is now a number.] Now, to get the correct c 2 , we sum 
over all possible intermediate states as we do for the case of an inclusive process and we get, 



x 2 (l-x) 2 p 2 2(27r) 3 fc'+- 



where, 



Here we have used, 



a.0 



2(27r) 3 fc'+ 
P = 



}— g 2 r(p,p',k')T^{p,p',k') 



N 2 -l 
2N 



1 



]p+][l-x] { [x] 
1 



\p+][l-x] {[x] 



0h % 



N 2 -l 



Ef rjiCl rpa \ 

aft 

E X\2X\ 2 = i 

E'A'l 



(3.9) 

(3.10) 
(3.11) 

(3.12) 
(3.13) 
(3.14) 



Now, 



pipif — 



=> c^{x, k±) - 

as k' + = xp. Integrating over k^, wc get, 



p 2 (l — x) 2 I X 2 X 
[l + (l-x) 2 ] 



2k 2 ± 



x 2 (l — x) 2 p 2 



"1 „2 



16n 3 p" 



N 2 -l\ l + (l-.x) 2 2 



2N 



(3.15) 
(3.16) 
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" ' 16ttV V 2N 

c? a „ , /A x 2 



where, 



and a = A and /i are some ultraviolet and infrared cutoff respectively. 
P qq can be found using the relation 



'jSf 2 -i\ l + x 2 

/>„„<.<> :) = (__)_Z_. (3.19) 



B. Calculation of P q c and Pgg 

Now, if we want to calculate P q c or Pgg, we have to start with a gluon state which is given by, 



\gluon) = \ip{k)) 



2„± r 



+ E / / d j==^ P - « - A 2 )|0) C2 („ *) 

1 W /" g3 (fc _ fci _ fc 2 )qt( fcl , Axjqt^, A 2 )10)4(fc 1; fc 2 ) (3.20) 

Then we proceed along the same way as before. Consider the eigenvalue equation for the gluon state and take the 
projection onto a quark- antiquark state and two-gluon state separately. 

\quark - antiquark) = \p'q') = 1 hUp' >! X )S (q 1 , A' 2 )|0) (3.21) 

v 2 ( 2?r ) 3 

\two - gluon) = \k[, k' 2 ) = 1 1 a^{k[, Ai)a f (^, X' 2 )\0) (3.22) 

^2(2nfk+' ^2(2^+' 



In the first case, 1st and 2nd term from the gluon state will contribute and we get as in equation (3. 

ci-^S^ = *vip^>i 2 = ^Sn^™ 2 - ^ 

Here, we have used the following momentum parametrization. 

total gluon mom. k = (k, k, 0j_) 

/ k 2 \ 

antiquark mom. q' = I (1 — x)k H — - 1 , (1 — x)k, —kj_ I 

\ 2(1 — x)k J 

( k 2 \ 

gluon mom. p' = i/cH — , xk, k± ) . (3.24) 

* \ 2xk J 

Now, with this parametrization (p'q'\P~\k) is exactly identical to (p'q'\P~\p) we already have calculated, with the 
difference that P now becomes 



G 



Now, summing over all possible intermediate states, we get, 

E \(p'a'\P v -\k)\ 2 = g 2 E (Wa) E (xlr'xx.xl^xx,) $ J 



A C A 



/3a AiA 2 

i*J 
A 6 A 



= 5 9 Q) 21*1^*4 (3.27) 

[The factor 2 comes from the summation over Ai. ] 

At this point we assume that the result is independent of initial gluon polarizations and hence we add the result 
for different polarization and divide by two to get the final result. 

-El(PV|^|fe)| 2 = ^Wt (3 . 28) 



Now, 



k 2 (l — x) 2 \ x x 2 J k 2 x 2 (l — x) 



cl{X ^ ] ^i&k\ 92[x2 + [l ^ X)2] ^i (3 ' 30) 



Integration over k± gives, 



where, 



c 



P qG = l[x 2 + (l-x) 2 ]. (3.32) 



In the second case, when we take t he p rojection onto two-gluon state, 1st and 3rd term in equation (3.20) will 



contribute, and we get as in equation (3.- 



x 2 (l — x) 2 k 2 



2 



2(2vr) 3 k 3 {l -x)x 



\(k[k' 2 \p-\k)\ 2 . (3.33) 



[Again, we have extracted the factors arising from the normalization of states and the energy-momentum conserving 
(5-function, and write them separately so that the remaining factor is just a number.] 

Now, summing over all the intermediate gluon states and averaging over the initial gluon states, we get the following 
result. 



J2\(q,k\Pv\p)\ 2 = 9 2 £/ o6 7.fc E riilTti ' j ' 1 ' 



be 



' i j I 
i> f V 

^EC4E«E^< (3-34) 

Ai A2 A3 
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Now using the relations 



be 



EJ* J' — X.., 

A 2 



J2\(q,k\P-\ P )\ 2 = g 2 Nlj2 TijlT 



Now, 



1 1 

.x 2 (1 — x) 2 



J2\(i,k\p-\p)\ 2 =g 2 N 



i i 



x 2 (l-x) 2 



Thus, using this result, we get, 



x(l - x) + + 



x 1 — X 



Integration over k± gives, 



where, 



1 — x x 



■a c i ~ a , f A 



Pgg - 2Ar 



ar(l - a;) + 



x 1 — x 



(1 — x) x 



(3.35) 
(3.36) 

(3.37) 

(3.38) 
(3.39) 

(3.40) 

(3.41) 

(3.42) 



IV. CONCLUSION 



Using the normalization condition we can write c 2 = 1, in the equations ( |3.17 ), ( 3.31 ), and ( |3.41 ), since we are 
working only with the lowest order contributions in the coupling constant. Thus, we can see directly tha t the sp littin g 
f unct ions ar e pro portional to the corresponding probabilities as is the case in AP-equation. Equation ( 3.18 ), ( 3.19 ), 
(3.32) and (3.42) show the expressions for the splitting functions which are exactly identical to that derived in AP- 



paper. 

Of course, we have not considered the contributions coming from the infrared singularity (x — > 1 limit) in the 
expression of splitting function. One usually puts them on the basis of some physical arguments or by direct loop 
calculations. 

In principle one should be able starting from the first principle, to derive the AP-cquation in the hamiltonian 
framework. The very success of our calculation in obtaining the splitting functions suggests that it may be plausible 
to try and extrapolate it to obtain AP-equation, which is, of course, a non-trivial job. 
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